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Abstract — The temperatures of prolate and oblate spheroidal dust grains in the envelopes 
of stars of various spectral types are calculated. Homogeneous particles with aspect ratios 
a/b < 10 composed of amorphous carbon, iron, dirty ice, various silicates, and other materials 
are considered. The temperatures of spherical and spheroidal particles were found to vary 
similarly with particle size, distance to the star, and stellar temperature. The temperature 
ratio Td (spheroid) /Td (sphere) depends most strongly on the grain chemical composition and 
shape. Spheroidal grains are generally colder than spherical particles of the same volume; 
only iron spheroids can be slightly hotter than iron spheres. At a/b ~ 2, the temperature 
^ ■ differences do not exceed 10%. If a/b ^ 4, the temperatures can differ by 30-40%. For a 

*p"( fixed dust mass in the medium, the fluxes at wavelengths A <; 100 /im are higher if the grains 

are nonspherical, which gives over estimated dust masses from millimeter observations. The 
effect of grain shape should also be taken into account when modeling Galactic-dust emission 
properties, which are calculated when searching for fluctuations of the cosmic microwave 
background radiation in its Wien wing. 

•»— i . 
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1. INTRODUCTION 

The observed infrared and submillimeter emission from interstellar clouds, circumstellar en- 
velopes, and galaxies is generally thermal emission of dust heated by stellar radiation or shock 
waves. When infrared spectra of these objects are computed, the dust temperature must be cal- 
culated. This temperature is also used to determine the mass and thermal balance of the matter 
in various objects and is important for the formation of molecules on the grain surfaces. 

Calculations of the interstellar dust temperature were initiated in the 1940s [see van de Hulst 
(1949) for a discussion]. The equilibrium temperature of spherical dust grains is commonly consid- 
ered (see, e.g., Mathis et al. 1983). However, it has been known for fifty years [since the discovery 
of interstellar polarization by Hiltner (1949), Hall (1949), and Dombrovskii (1949)] that there are 
nonspherical aligned particles in the interstellar medium. 



1 To the memory of Gennadii Borisovich Sholomitskii, an enthusiast for research in the field of infrared and 
submillimeter astronomy. 
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Nonspherical particles appear to be also present in the circumstellar dust shells. The variations 
in the position angle of linear polarization with time and wavelength observed in red giants provide 
circumstantial evidence for this (Dyck and Jennings 1971; Shawl 1975). Even after correction for 
the interstellar polarization, the position-angle difference in the blue and in the red can reach 
20° — 60°. This behavior is very difficult to explain in terms of the model of a single star with 
a shell containing spherical particles alone. If, alternatively, there are nonspherical grains in the 
shell, then variations in the degree and direction of grain alignment may result in the observable 
variations of the polarization angle. 

The first attempt to take into account the effect of the shape of interstellar grains on their 
temperature was made by Greenberg and Shah (1971). These authors considered metallic and 
dielectric Rayleigh spheroids and infinite ice cylinders of 0.1-/im radius. They concluded that 
nonspherical particles were approximately 10 % colder than spheres, a result that entered the 
books on interstellar dust (Whittet 1992). 

Recently, Fogel and Leung (1998) have computed the infrared radiation of fractal dust grains 
produced by two processes of stochastic growth and composed of amorphous carbon and silicate. 
They concluded that the temperature of nonspherical particles was typically 10-20 % lower than 
that of spherical ones, which results in a longward shift of the maximum of the object's radiation. 
Fogel and Leung (1998) proposed to consider the fractal particle size as a parameter of the grain 
shape, but they ignored the effects of alignment. 

Here, we study in detail the dependence of the temperature of spheroidal circumstellar dust 
grains on their shape and alignment. We consider prolate and oblate particles of various sizes, 
which are composed of a variety of absorbing and dielectric materials and which lie at various 
distances from stars with various temperatures. The nonsphericity effect of interstellar dust grains 
on their temperature was discussed by Voshchinnikov et al. (1999). 

2. THE MODEL 

2.1. The Radiation Field 

Dust particles are present in the envelopes of late-type stars (red giants and supergiants) and 
hotter stars, such as Herbig Ae/Be stars. Dust grains are heated mainly by the absorption of 
stellar radiation. As was pointed out by Lamy and Perrin (1997), in some cases, it is important 
to take into account the star's true spectral energy distribution as well. In addition, in the outer 
regions of optically thick dust shells, the maximum in the energy distribution is redshifted (see, 
e.g., Bagnulo et al. 1995). For simplicity, however, we assume the energy distribution to be a 
blackbody one with an effective temperature T*. Since we are going to compare the temperatures 
of particles under the same conditions, this assumption is of no fundamental importance. 

In most cases, we take the stellar temperature to be T* = 2500 K, typical of late-type giants 
and supergiants (Pegourie 1987; Lorenz-Martins and Lefevre 1995). The effects of variations in 
T* are discussed in subsection 4.5. 

2.2. Dust Particles 

Chemical composition. Particles of amorphous carbon and amorphous silicates are most 
commonly considered as the major sources of the infrared radiation observed from carbon and 
oxygen stars, respectively. Evidence in support of these materials follows both from theoretical 
models of dust formation (Gail and Sedlmayr 1984) and from laboratory experiments (Jager et 
al. 1994). 



The specific type of silicate or carbon material in the circumstellar medium is very difficult 
to determine. In addition, there is most likely a mixture of various dust components in the 
shells simultaneously. For example, bands of silicon carbide, sulfide silicates, and even amorphous 
silicates were detected in the spectra of some carbon stars (Baron et al. 1987; Goebel and Moseley; 
Little-Marenin 1986). Several emission bands found in the spectra of carbon stars were identified 
with crystalline silicates (Waters et al. 1999). Finally, iron and oxide particles can apparently 
condense in circumstellar envelopes irrespective of the C/O ratio. 

In our modeling, we used the six materials that were previously chosen by Il'in and Voshchin- 
nikov (1998) when considering the effect of radiation pressure on dust grains in the envelopes of 
late-type stars: amorphous carbon, iron, and magnetite (FesC^) as examples of strongly absorb- 
ing materials, as well as astronomical silicate (astrosil), transparent glassy pyroxene, and artificial 
dirty silicate (Ossenkopf et al. 1992; OHM silicate) as silicates of various types. References to 
the papers from which we took the optical constants of these materials can be found in Il'in and 
Voshchinnikov (1998).Q This set of materials was extended to include carbon material (cellulose), 

o 

which was produced by pyrolysis at the temperature of 1000 C (cellOOO; Jager et al. 1998), and 
dirty ice, which was used in the classical study by Greenberg and Shah (1971). In the latter case, 
the imaginary part of the refractive index was chosen to be k = 0.02 in the wavelength range 
0.17-1.2 /im, as was done by Greenberg (1970, 1971). 

Shape. The formation of only spherical particles in the envelopes of late-type stars has been 
considered thus far (see, e.g., Draine 1981; Gail and Sedlmayr 1985; Fadeyev 1987; Fleischer et al. 
1992; Cadwell et al. 1994). The theory of nucleation and growth of nonspherical particles is still 
at the initial stage of its development. 

We assume the circumstellar dust grains to be prolate and oblate homogeneous spheroids with 
aspect ratios a/b (a and b are the spheroid semimajor and semiminor axes, respectively). By 
varying a/b, we can model the particle shape over a wide range: from spheres to needles and 
disks. 

Size. Dust grains form and grow in the envelopes of late-type stars. They range in size 
from tiny particles to particles with radii up to 1 /mi or more [see Lafon and Berruyer (1991) 
for a discussion]. The upper limit of the grain size distribution is uncertain and is the subject of 
debate. However, it follows from model calculations that, in general, the particle size in oxygen 
stars is larger than that in carbon ones (Jura 1994, 1996; Bagnulo et al. 1995). 

In order to compare the optical properties of particles of the same volume but different shape, 
it is convenient to characterize the particle size by the radius ry of a sphere equal in volume to a 
spheroid. The spheroid semimajor axis is related to ry by 



a 
,V l b 



x 2/3 

) « 



for prolate spheroids and by 



a = r v { l ) (2) 

for oblate spheroids. In our calculations, we considered particles with ry = 0.005 — 0.5 /im. 

Structure. The dust grains growing in circumstellar shells can be fluffy or porous. To model 
the effect of porosity, we used Bruggeman's rule (Bohren and Huffman 1986) and obtained the 
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Data on the refractive indices can also be extracted from an electronic database of optical constants (Henning 



et al. 1999) via Internet at http://www.astro.spbu.ru/JPDOC/cntry.html.) 



mean effective dielectric function e c g of an aggregate composed of n materials with dielectric 
functions £$, 

£?; — £eff 



£/, 



o. 



(3) 



i=i £ « + 2^cff 

where /j is the volume fraction occupied by the material of type i. The temperature is calculated 
for compact particles with e e s. We considered spheroids composed of vacuum (e = 1) and cellulose 
with vacuum fractions from to 0.9. 

Orientation. Collisions of dust grains with atoms and molecules cause rapid grain rotation 
with angular velocities > 10 5 s _1 . Interstellar particles are believed to rotate around the direction 
of maximum moment of inertia, and, in general, the angular momentum is parallel to the magnetic 
field (Spitzer 1981). Circumstellar particles can be aligned by anisotropic radiation or gas fluxes 
(Dolginov et al. 1979). Radial grain motion in the shells must apparently cause particle rotation in 
the planes containing the radius vector. However, nonradial gas flows or the helical circumstellar 
magnetic fields produced by stellar rotation (see, e.g., Woitke et al. 1993) can also result in a 
different grain alignment. 

In our modeling, we considered two types of grain orientation: particles randomly oriented 
in space (3D orientation) and in a plane (2D orientation or complete rotational orientation). In 
the latter case, the major axis of a rotating spheroid always lies in the same plane. The angle 
Q between the particle angular velocity and the wave vector of the incident radiation is a model 
parameter (0° < fi < 90°). 

3. BASIC EQUATIONS 

Let us consider a dust grain at distance R from a star of radius R* and temperature T*. 
The stellar radiation is assumed to be unpolarized. The equilibrium grain temperature T& can 
be determined from Kirchoff 's law by solving the energy balance equation for the absorbed and 
emitted energy (erg s _1 ) 



W / C abs (A) 7r5 A (T*) dA = / C em (A) nB x (T d ) dA , (4) 

JO Jo 

where C a t, s (A) and C cm (A) are the orientation-averaged absorption and emission cross sections, 
-kB\(T) is the blackbody flux with temperature T (erg cm~ 2 s _1 /im" 1 ), and W = R\jR 2 is the 
radiation dilution factor. 

For particles randomly oriented in space, the absorption cross sections must be averaged over 
all orientations: 

Cabs = / o < 5abs I ( m A,rv,A,a/6,a) + Qj b E s (m A ,r v ,A,a/6,a) G(a)smada. (5) 
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Here, m\ is the refractive index of the grain material [rn\ = e x ), a is the angle between the 
spheroid rotation axis and the wave vector (0° < a < 90°), and G is the geometric cross section 
of the spheroid (the area of the particle shadow): 
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for prolate spheroids and 
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for oblate spheroids. 

In the case of complete rotational orientation, the absorption cross sections are averaged over 
all rotation angles <p. For prolate spheroids, this yields 



*-<"> = If 5 
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<2abs ( m A, TV, A, a/b, a) + Q abs (m A , r v , A, a/b, a) G(a)d(j) 



where the angle a is related to Q and <fi (cos a = sin Vt cos <ft). For oblate spheroids arbitrarily 
oriented in a plane, we have Q = a and 



-^ 2D, 



C£(n) = 2 [Q™("»A,rv, A,a/6,n) + Q™(m x ,r v , \,a/b,Q)\ G(Q). 



(9) 



The energy emitted by a particle is proportional to its surface area. The emission cross section 
can then be calculated as follows: 
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where 
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for prolate spheroids and 



S = 2rrrl 
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(10) 
(11) 
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for oblate spheroids, e = J 1 — (a/6) 2 . 

In Eqs. (5), (8), (9), and (10), the superscripts TM and TE refer to two polarizations of the 
incident radiation (TM and TE modes). The efficiency factors Q a b s ' can be calculated using an 
exact solution to the problem of diffraction of a plane electromagnetic wave by a homogeneous 
spheroid by separation of variables method [see Voshchinnikov and Farafonov (1993) for more 
detail). We used the benchmark results of calculations of the efficiency factors from Voshchinnikov 
et al. (2000) to thoroughly test the computational program. 

4. RESULTS AND DISCUSSION 

In this section, we discuss the dependence of the spheroidal-grain temperature on model pa- 
rameters: the particle chemical composition, size, shape, and orientation, as well as the stellar 
temperature and distance to the star. We also consider the effects of porosity and partial polariza- 
tion of the incident radiation on the particle temperature. Since the main goal of our study is to 
analyze the effects of change in the particle shape, we compare the temperature ratio of spheroidal 
and spherical particles of the same volume. 

4-1. A Graphical Method of Temperature Determination and a Blackbody Model 

This method, proposed by Greenberg (1970, 1971), allows the causes of temperature differences 
between dust grains with different characteristics to be clearly established if the particles are 
assumed to be in an isotropic radiation field. Their temperature can then be determined by 
solving the following equation of thermal balance: 



W / C abs (A)4vr J B A (T,)dA 
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Figure 1: The power absorbed (emitted) by dust particles with ry = 0.01 /im in an isotropic 
radiation field. The straight line corresponds to a blackbody. The curves refer to spherical and 
prolate spheroidal (a/b = 10, 3D orientation) cellulose (cellOOO) particles. A graphical method 
of determining the grain temperature is shown for T* = 2500 K and W= 10~ 8 . In this case, the 
temperatures of spherical and spheroidal particles and a blackbody are respectively, Td (sphere) = 



76.8 K, T d 3D 



48.7 K, and T d BB 



2500(KT 8 ) 1 / 4 = 25 K. 



If the spheroidal particles are randomly oriented in space (3D orientation), then numerical esti- 
mates show that 



■p=y3D ^ L/ em 
^abs ~ ~~T~ 



(14) 



The integrals on the left- and right-hand sides of Eq. (13) then depend on the temperature T 
alone (for given particle chemical composition, size, and shape). 

Figure 1 shows the power for the absorbed (emitted) radiation (in erg s _1 ) for spheres and 
spheroids of the same volume (ry = 0.01 /im) composed of cellulose (cellOOO). The method of 
temperature determination is indicated by arrows: from the point on the curve corresponding to 
the stellar temperature (T* = 2500 K), we drop a perpendicular whose length is determined by 
the radiation dilution factor and then find the point of intersection of the horizontal segment with 
the same curve for the power. This point gives the dust grain temperature T d determined from 
the equation of thermal balance (13). 

It follows from Fig. 1 that the different temperatures of spheres and spheroids result from 
different particle emissivities at low T. The differences are most noticeable for cellulose. For other 
materials, the pattern of dependence is preserved: the emissivity of spheroidal particles at low 
temperatures is larger than that for spheres and closer to the blackbody one (see Fig. 2). Only 
iron particles constitute an exception: spheres composed of them at T ^ 30 K emit more energy 
than spheroids with a/b ~ 2 — 6 but less energy than spheroids with a/b ^ 8 (Fig. 3). 

The straight lines in the upper parts of Figs. 2 and 3 refer to a blackbody, for which the 
emissivity is proportional to T 4 , and the absorption efficiency factors are equal to unity, 
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Using Eq. (15) and Stefan-Boltzmann's law, instead of Eq. (13) we obtain for the 3D orientation 

(16) 
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Figure 2: The power absorbed by dust particles with ry = 0.01 /im in an anisotropic radiation 
field. The curves refer to spherical, prolate, and oblate spheroidal amorphous carbon (AC1) 
particles with a/b = 10. The straight line corresponds to a blackbody. 
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Figure 3: The power absorbed by iron dust particles with ry = 0.01 /im in an anisotropic radiation 
field. 



Considering that the mean projection area of any convex figure is a quarter of its surface area 
(see, e.g., Hildebrand 1983) fl 
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4 ' 
we have for spheres and spheroids randomly oriented in space 

T d BB = T,iy 1 / 4 . 

For the 2D orientation, the blackbody temperature is 
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Expressions for the mean cross sections in the case of complete rotational orientation were derived 
by Il'in and Voshchinnikov (1998): 
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for prolate spheroids and 



G 2D (fi) = G(Q) 



(20) 



(21) 



for oblate spheroids. In Eq. (20), the complete elliptic integral of the second kind is denoted by 
E{m). Based on Eqs. (20) and (21), we can easily determine the blackbody-temperature ratios 
for the two extreme cases of 2D orientation (fi = 0° Q = 90°): 
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T BB (fi = 90°) ~ Kb J (23) 

for oblate spheroids. We can also calculate the temperature ratio of a spheroid for the 2D orien 
tation and a sphere (or a spheroid for the 3D orientation): 



T d BB (fi) 



T d BB (sphere, 3D) 
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(24) 



Note that Eqs. (22)-(24) allow the geometric effects on the temperature of nonspherical dust 
grains to be estimated. For example, T BB (fi = 0°)/T BB (sphere) = 1.06 (1.18) for a/b = 10 and 
T d BB (fi = 90°) /T d BB (sphere) = 0.95 (0.66) for prolate (oblate) spheroids. 

4-2. Dependence on the Distance from the Star 

The temperature of particles of any shape decreases with distance from the heating source 
(star). However, the temperature ratio of nonspherical and spherical particles varies over a narrow 

3 Eq. (17) for prolate and oblate spheroids can be proved by directly integrating Eqs. (6) and (7). 




0.85 
1.05 



-d 
o 



si 



prolate, AC1 



vvvvv 2D, n = 

AAAAA 2D, = 90° 

®m ®®m 3D 



oblate, AC1 



-SJ ^- 




10 



100 



1000 10000 100000 

R/R* 



Figure 4: The temperature ratio of prolate and oblate spheroidal (a/b = 4) and spherical 
amorphous carbon (AC1) grains with ry = 0.01 /im at various distances from a star with T* = 
2500 K. The sphere temperature is T d (sphere) = 780.6 K, 318.0 K, 130.8 K, 52.3 K and 20.8 K 
for R = 10, 10 2 , 10 3 , 10 4 and 10 5 R*, respectively. 



range. As follows from Fig. 4, which shows the results of our calculations for amorphous carbon 
particles, this ratio decreases approximately by 5 % as the distance to the star increases from 
10 to 10 5 R±. In this case, the sphere temperature decreases from 780.6 to 20.8 K. Given the 
insignificant variation of Td(spheroid)/Td(sphere) with R, below we consider the results only for 
R=10 4 R*(W = 10~ S ). 

4-3. Dependence on the Grain Size 

The temperature variations of carbon and silicate dust grains with particle radius are given 
in Tables 1 and 2. As ry increases, the particle temperature slightly rises and then begins to 
fall, as was also noted previously (Greenberg 1970). The temperature of spheroidal particles 
is always slightly lower than that of spherical ones; the difference can reach 10 % at a/b = 4. 
The temperature difference between spheres and spheroids is largest for the 2D orientation and 
Q = 90°. Note also that, since T& for spheres and spheroids vary with ry similarly, the temperature 
ratio is almost independent of ry. 

4-4- Dependence on the Grain Shape and Chemical Composition 

The effects of the chemical composition and shape of dust grains on their temperature are most 
significant compared to the variations of other model parameters. These effects are illustrated in 
Figs. 5 and 6, which show the temperature ratios of spheroids and spheres for particles with 
metallic and dielectric properties, respectively. The variations of T^ (spheroid) /T& (sphere) have 
the following general tendency: this ratio decreases with increasing a/b and absorptive material 
properties, with the effect increasing. Besides, the following inequality holds for particles of the 
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Table 1: The temperatures (in K) of spherical and spheroidal (a/6 
T* = 2500 K, R = 10 000 R* 



4) amorphous carbon grains; 
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Table 2: The temperatures (in K) of spherical and spheroidal (a/b = 4) astronomical silicate 
grains T* = 2500 K, R = 10 000 R* 
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Figure 5: The temperature ratio of prolate and oblate spheroidal and spherical grains composed 
of various materials with metallic properties; T* = 2500 K, ry = 0.01 /an, and R = 10 4 R*. The 
sphere temperature is Td (sphere) = 52.3 K, 63.0 K, 59.2 K, and 120.0 K for particles of amorphous 
carbon, magnetite, cellulose, and iron, respectively. 
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(25) 



The difference between prolate and oblate particles, which is noticeable in Figs. 5 and 6, is 
determined by different optical properties of the particles [see Voshchinnikov and Farafonov (1993) 
for more detail]. If we restrict ourselves to a/b ^ 4, then it is easy to see that the dielectric and 
metallic spheroids are, respectively, 10 % and 20 % colder than the corresponding spheres. The 
temperatures of prolate and oblate particles can differ by 30-40 % from the temperature of spheres 
(Fig. 5). If we further increase the aspect ratio a/b, then the temperature gradually ceases to 
drop. For example, for prolate cellOOO particles, T d BD /T d (sphere) = 0.61, 0.53, 0.50, and 0.50 for 
a/b = 10, 20, 50, and 100, respectively 

The temperature behavior for iron particles is peculiar: for a/b ^ 8, the spheroids are 10-20 % 
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Figure 6: The same as Fig. 5 for particles with dielectric properties. The sphere temperature is 
Td(sphere) = 34.9 K, 24.0 K, 42.4 K, and 33.7 K for particles of astrosil, pyroxene, OHM silicate, 
and dirty ice, respectively. 
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Figure 7: The temperature ratio of prolate and oblate spheroidal (a/b = 4) and spherical 
amorphous carbon (AC1) grains with ry = 0.01 /zm at R = 10 4 R± from the star. The sphere 
temperature is T d (sphere) = 10.5 K, 20.4 K, 63.1 K, 106.4 K, and 226.4 K for T* = 500 K, 1000 K, 
3000 K, 5000 K, 10 000 K, respectively. 



hotter than the spheres. The reason for this is seen from the behavior of the particle emissivity 
shown in Fig. 3. 

4-5. Dependence of the Stellar Temperature 

Nonspherical particles may be present in the shells of stars of various spectral types and in 
the interplanetary medium. The effects of variations in the stellar temperature are shown in Fig. 
7 for amorphous carbon particles with ry = 0.01 /im. Low effective temperatures (up to 500 K) 
correspond to possible conditions in the outer parts of optically thick circumstellar shells. It follows 
from Fig. 7 that there is no effect for T,± (spheroid) /Td(sphere), although the absolute temperature 
varies over a fairly wide range: from Td(sphere) = 10.6 K at T* = 500 K to Td(sphere) = 226.4 K 
at T* = 10 000 K. 

4-6. Porous Dust Grains and Polarized Incident Radiation 

An increase in the fraction / of vacuum in porous particles causes a decrease in the effective 
refractive index, i.e., a reduction in particle absorptivity. In this case, the shape effects gradually 
disappear, and the temperature of spheroidal particles approaches the temperature of spheres. 
This effect is illustrated in Table 3, which shows that Td(spheroid)/Td(sphere) increases consider- 
ably for cellulose particles at / ^ 0.5. 

The aligned nonspherical particles in the inner parts of circumstellar dust shells can polarize the 
stellar radiation. In principle, the local linear polarization of scattered radiation can be significant. 
The absorption efficiency factors for the TM and TE modes must then enter into the absorption 
cross sections [Eqs. (5), (8), and (9)] with different weights, depending on the polarization P of the 
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Table 3: The temperatures (in K) of porous spherical and prolate spheroidal (a/b = 10) cellulose 
(cellOOO) grains; r v = 0.01 , T* = 2500 K, R = 10 000 i?* 



1 — / T d (spheroid) /T d (sphere) T d (sphere) 





2D, tt = 0° 


2D 


, tt = 


90° 


3D 




1.0 


0.65 




0.59 




0.61 


59.2 


0.9 


0.65 




0.59 




0.61 


57.7 


0.7 


0.68 




0.63 




0.65 


53.5 


0.5 


0.77 




0.73 




0.74 


47.3 


0.3 


0.94 




0.91 




0.92 


45.0 


0.2 


0.98 




0.96 




0.96 


49.4 


0.1 


0.98 




0.97 




0.97 


54.8 



incident radiation. This causes an increase in the temperature of prolate spheroids if the electric 
vector of the incident radiation is parallel to the particle major axis (in this case, Q™ > Q]^s) and 
a decrease in their temperature otherwise (Q™ < Q™)- The reverse is true for oblate spheroids. 
However, the effect is not too large: even at P = 50 %, the temperature of prolate spheroidal 
(a/b = 10) cellulose particles varies between 33.2 and 38.3 K (T d D = 36.0 K at P — %). 

The following fairly general conclusion can be drawn from what was considered above: the 
temperatures of spherical and nonspherical dust grains are proportional, with the proportionality 
coefficient being approximately the same for particles of different sizes at different distances in 
the envelopes of stars of diverse spectral types. This coefficient is determined only by the particle 
shape and by absorptive properties of the material of which they are composed. 

5. ASTROPHYSICAL IMPLICATIONS 

The dust temperature is an important parameter that determines the infrared spectra of various 
objects. The efficiency of grain growth and destruction, formation of molecules on the grain 
surfaces, and the alignment of nonspherical particles depend on T d (Voshchinnikov 1986). Besides, 
T d determines the polarized submillimeter emission of nonspherical particles (Onaka 2000) and 
the gas cooling rate in very dense interstellar clouds (Whitworth et al. 1998). 

The infrared flux at wavelength A emerging from an optically thin medium is proportional to 
the total number N of dust grains in the medium, the Planck function, which depends on the 
particle temperature T d , and the emission cross section C era (\) : 

Fm(X) = N ^j^-B x (T d ) , (26) 

where D is the distance to the object. The quantities T d and C em (A) in Eq. (26) depend on the 
particle shape. 

If we assume the chemical composition and sizes of all particles in the medium to be the same 
and if we change only the grain shape, then, when spheres are replaced by nonspherical particles, 
the position of the maximum in the spectrum of thermal radiation is shifted longward, because 
the temperature of the latter is higher (Fig. 8). 
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Figure 8: The normalized fluxes emerging from a medium containing the same (by mass) amount 
of spherical and prolate spheroidal (a/b = 10, 3D orientation) cellulose (cellOOO) particles; ry = 
0.01 /zm, T* = 2500 K, and R = 10 4 R±. The particle temperatures are Td(sphere) = 59.2 K and 
T d 3D = 36.0 K. 
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Figure 9: The same as Fig. 8 for the flux ratio. The ratios of the emission cross sections and 
fluxes for a blackbody are also shown. The latter are determined by the particle temperature 
alone. 



The flux ratio is specified as follows: 

^spheroid / \ \ 



IR 

j-isphcrc 



(A) 



C: P m hcroid (A) ff A [T d (spheroid)] 
CX hcrc (A) 



B\[T A (sphere) 



(27) 



It is shown in Fig. 9, where the ratios of the cross sections and the Planck functions in Eq. (27) 
are also displayed. Despite the different particle temperatures, the Planck functions differ insignif- 
icantly. The main differences in the fluxes result from different emission cross sections C em (A): 
they are considerably larger for nonspherical particles at A > A max (Fig. 9). This effect was 
previously noted by Ossenkopf and Henning (1994). 

If the dust mass in an object is determined from the observed millimeter flux, then the 
Rayleigh- Jeans approximation can be used for the Planck function. For the same flux Ftr(A), the 
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dust mass depends on the particle emission cross sections and temperatures, while the ratio 

Mf here C^ lcroid (A) T d (spheroid) 

M sphcroid " C - e Sphere (A) Td (s P here) 



(28) 



shows the extent to which the values of M& differ when changing the grain shape. 

The dust mass in galaxies and molecular clouds is commonly estimated from 1.3- mm observa- 
tions (Siebenmorgen et al. 1999). At this wavelength, the ratio of the cross sections for cellulose 
particles with a/b = 10 (see Figs. 8 and 9) is ~ 50, and the temperature ratio is 36.0/59.2 = 0.61, 
which gives approximately a factor of 30 larger dust mass if the particles are assumed to be 
spheres rather than spheroids. In other (not so extreme) cases, the dust-mass overestimate is 
much smaller. For example, for prolate amorphous carbon particles, M^ p ere /Mj P ero1 ss 1.2,2.3, 
and 6.4, for a/b = 2, 4, and 10, respectively. 

The degree of alignment of interstellar and circumstellar dust particles for rotational (Davies- 
Greenstein) orientation depends on the dust-to-gas temperature ratio: the smaller is this ratio, 
the more effective is the alignment of nonspherical particles (Greenberg 1970; Voshchinnikov 
1986). In this case, as follows from Tables 1 and 2 and from the discussion in subsection 4.4, the 
nonspherical particles have a minimum temperature when the magnetic field is perpendicular to 
the line of sight (Q = 90°), i.e., when the polarization of the forward transmitted radiation is at 
a maximum. However, this remark more likely refers to protostellar envelopes than to interstellar 
particles, where the dust grains are heated by isotropic radiation. 

Yet another important astrophysical process is the formation of molecules on the grain surfaces. 
Since the efficiency of this process strongly depends on the grain temperature, a lower temperature 
of nonspherical particles can significantly facilitate the formation of molecules on their surfaces 
[see Voshchinnikov et al. (1999) for a discussion]. 

Finally, it should be mentioned that the sublimation (along with nucleation) of spherical and 
nonspherical particles must take place at different distances from the star. If we consider the 
evaporation of ice cometary grains in the interplanetary medium, then, despite the weak shape 
effect for ice (see Fig. 6), the differences for spheres and spheroids prove to be noticeable enough. 
Ice spheres reach the sublimation temperature [T d (sphere) = T su bi = 100 K] and begin to evaporate 
at a distance R ~ 20 AU from the Sun. For prolate and oblate spheroids (a/b = 10, 3D orientation), 
this takes place closer to the Sun: at R ~ 18 and 16 AU, respectively. 

6. CONCLUSIONS 

Nonspherical interstellar and circumstellar dust grains must be slightly colder than spheres of 
the same volume composed of the same material. An analysis of the shape effects for spheroidal 
particles shows that the temperature differences do not exceed 10 % for a spheroid aspect ratio 
a/b rs 2. If a/b > 4, then the temperature differences reach 30-40 %. Only iron spheroids can be 
slightly hotter than spheres. All shape effects are mainly associated with different emissivities of 
spherical and nonspherical particles at low temperatures. Grain porosity causes a reduction in the 
shape effects, whereas polarization of the incident radiation can either reduce or enhance them. 

The calculated fluxes at A ^ 100 /zm turn out to be higher if the grains are assumed to be 
nonspherical, resulting in an overestimation of the dust mass inferred from millimeter observations. 

Finkbeiner and Schlegel (1999) point out that the Galactic-dust thermal emission at long 
wavelengths gives a major contribution to the observed background infrared radiation and must 
be subtracted in the search for fluctuations of the cosmic microwave background radiation in its 
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Wien wing. In this case, very small (in magnitude) effects compared to which the effect of grain 
shape is significant and must be taken into account in modeling are sought. 
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